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OPTIMAL LINK SCHEDULING AND POWER CONTROL 
FOR CDMA WIRELESS MULTI-HOP NETWORKS 

Consider a generic multi-hop network with a set of wireless devices (abstracted as 
nodes) communicating with each other in a certain pre-specified configuration (abstracted 
as directed links). Nodes could represent laptops, PDA's, Access Points etc., while 
directed links indicate data transfer between nodes. AH nodes share a common frequency 
spectrum and use DS-CDMA, transmitting their data using a signature sequence. This 
results in each link's transmission interfering with other ongoing transmissions in the 
network We have developed an algorithm that maximizes the weighted sum of long-term 
average link rates (throughput) of a multi-hop network when each link in the network is 
subject to peak and average power constraints. We call any subset of the links in the 
given configuration a transmission mode. Given the set of directed links and the signal 
attenuation between transmitters and receivers in the network, our proposed algorithm 
determines the throughput maximizing strategy. Specifically, it computes the optimal set 
of transmission modes, the relative duration for which each transmission mode is active 
and the power of signal transmission of every transmitter in each transmission mode. For 
more details refer to the attached manuscript. 

Using essentially the same methodology, an algorithm has been developed that 
imnimizes a weighted sum of the average power used by transmitters, subject to a lower 
bound on the average data rate on each link, and peak power constraints per transmitter. 
The algorithm yields a schedule and power allocations for each transmitter. Moreover, an 
additional output of the algorithm is a sensitivity value, for each link, which specifies the 
marginal increase in average (weighted) power required for a marginal increase in the 
average data rate on the given link. These link sensitivities can be used as link weights in 
a shortest path routing algorithm, enabling energy efficient routing. Although the 
methodology was conceived for DS-CDMA networks, the same methodology could be 
applied in wireless FDMA networks as well. 

Existing approaches pre-specify the link scheduling policy based on a graph- 
theoretic approach. They schedule links that are geographically far away (virtually 
interference-ftee) from each other simultaneously, e.g. maximum independent set of 
nodes in a graph. In our approach, we explicitly account for multi-access interference and 
de termine the throughput maximizing transmission strategy. In fact we find that as the 
level of ambient noise increases, the optimal policy schedules links in close spatial 
proximity simultaneously. A distinguishing feature of our approach is in the use of 
convex duality to solve this resource allocation problem. We establish that optimizing 
transmissions over the infinite time horizon can be achieved by optimizing over every 
single slot Next, we compute the optimal subsets of links by solving the reduced dual 
problem over a single slot We have developed an iterative descent algorithm that solves 
the dual problem in a finite number of steps. Solving the dual problem yields us a small 



set of transmission modes. Finally, we determine the optimal time-sharing of these 
transmission modes to maximize network capacity. 

Given a set of wireless devices and a certain pre-specified configuration of links, 
sharing a common frequency spectrum, our algorithm specifies a transmission strategy 
maximizes the weighted sum of long-term link rates when each link is subject to peak 
and average power constraints. The main input to our algorithm is the signal attenuation 
between all transmitters and receivers in the network. 

We solve the above constrained optimization problem using the theory of convex 
duality First, we establish that finding the best transmission strategy over the infinite 
time horizon can be achieved by optimizing over every single slot. Next, we compete the 
optimal transmission modes by solving the reduced dual problem over a single slot. We 
have developed an iterative descent algorithm that solves the reduced dual problem in a 
finite number of steps. Solving the dual problem yields us a small set of optimal 
transmission modes. Finally, we determine the optimal time-sharing of these transmission 
modes to maximize network capacity through a series of matrix inversions, each matrix 
sparse and comprising of either l's or O's. 
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optimal Link Scheduling and Power Control in CDMA 

Multihop Wireless Networks 



Rene L. Cruz & Arvind V. San than am 



Abstract — 

In this paper, we study the problem of radio resource al- 
location in wireless multihop networks. The main objective 
is to develop integrated link scheduling and power control 
policies to maximize the network throughput Transmitters 
in the network are subject to peak and average power con- 
straints. Subject to these constraints, we find optimal sub- 
sets of concurrent transmitters and associated power lev- 
els to maximize network throughput, taking into account 
multi-access interference. We find that when ambient noise 
levels are sufficiently high, it is generally desirable to en- 
able simultaneous transmissions, even if they are in close 
geographic proximity. We compare our results against pure 
CDMA to illustrate the significant gains of our approach. 
The results in this paper are applicable to multihop and sin- 
gle hop networks. 

Index Terms — CDMA, Multihop Networks, Resource Al- 
location. 



I. INTRODUCTION 

In recent years, there has been a great deal of interest 
and development in supporting data applications in multi- 
access wireless networks. Wireless devices are however, 
severely limited in bandwidth and energy and must use 
their resources judiciously. Resource allocation policies 
in currently deployed systems like IEEE802.il though 
simple and distributed in nature, perform poorly even un- 
der moderate loads. Centralized scheduling of node trans- 
missions can greatly improve network capacity, e.g see 
[3] [4]. Consistent with this, the IEEE802.il standard 
includes a point coordination function (PCF) which cen- 
trally schedules users. The High Data Rate (HDR) stan- 
dard also schedules transmissions centrally on the down- 
link channel. The benefits of centralized scheduling for 
multihop networks are equally rewarding if not more. In 
this paper, we extend link scheduling from one-hop net- 
works to the multihop paradigm. 

Hie capacity of a wireless channel is a function of the 
signal to interference plus noise ratio (SIR) of that chan- 
nel. For instance, the Shannon capacity of a link % over a 
frequency bandwidth W is C(i) = Wlog 2 (l + SIR(i)) 

The authors are with the Dept. of Electrical and Computer Engg., 
UCSD, USA. E-mail: asanthan@ucsd.edu . 



(assuming the noise is Gaussian). The dependence of link 
rates on powers naturally raises the possibility of improv- 
ing capacity through a suitable choice of powers. Before 
delving into specifics, we define some terms commonly 
used in this paper. A directed link / e {£, j} in a multihop 
network of N nodes facilitates data transfer from node i 
to node j. We call any subset of L a potential links in a 
network a transmission mode. There are 2 La transmis- 
sion modes for an L a link network. A list of transmission 
modes and the relative time-fractions for which they are 
"used" comprises a link scheduling policy. For each trans- 
mission mode, we also need to specify the power used by 
active transmitters, which constitutes a power control pol- 
icy. 

The conventional approach to scheduling links in a 
multihop networks is to allow nodes far away from each 
other to transmit simultaneously. Such a link scheduling 
policy ensures essentially interference-free transmissions, 
thereby achieving high instantaneous individual link rates. 
Nodes in the vicinity 1 of each other have an edge in the 
channel contention graph Gch- Computing a maximum 
independent set [1] in the channel contention graph re- 
quires global topology information, entailing centralized 
processing. Good approximation algorithms that allow for 
distributed implementation would naturally be the best re- 
course. One such example is the minimum-degree greedy 
algorithm [2], a heuristic that chooses nodes with mini- 
mum degrees in subgraphs. 

In our approach, we explicitly model the effects of sig- 
nal interference, allowing concurrent transmissions in a 
close geographic proximity. By globally optimizing the 
transmission modes and the associated power levels, we 
can achieve much higher throughputs than previous ap- 
proaches. Specifically, we solve the following resource 
allocation problem: 

Problem: Find an optimal link scheduling and power 
control policy to maximize the network throughput, when 
transmitters are subject to peak and average power con- 
straints. 

The main input to our algorithm is the path-loss ma- 
trix of the entire network. This matrix consists of N 2 



1 TVi>icallv within a 2-hoo distance of each other. 
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gain v9!Rs t one between each pair of nodes. In solving 
this problem, we only specify the relative time-fractions 
for which each transmission mode is active, not the ex- 
act time when it is "on". Viewing the time-fractions as 
weights, transmission modes can be scheduled over slots, 
using algorithms like deficit round robin. 

II. System Model 
Consider a system with N nodes in a CDMA environ- 
ment A subset of L a transmission links (among the pos- 
sible N{N - 1) links) constitutes a network topology, and 
is assumed to be prespecified. We wish to maximize the 
weighted sum of throughputs on these links, as will be de- 
scribed below in equation (3). For a given link I = {i, j}, 
the transmitter node i uses a signal power P(l). The path 
gain from node i to node j is given by G(i, j). The path 
gain G{i y j) models the effects of signal attenuation due to 
distance, channel fading and shadowing. We assume that 
the pith gains G(i, j) are constant. Denoting the trans- 
mitting and receiving nodes of link I by T(l) and R (I) re- 
spectively, the signal to interference and noise ratio (SIR 
or 7(0) for link / is given by, 



The average power consumed by each link Z, P avg (I) is 
the time-average over an infinite horizon : 

l m 

P<"*(i) = limsup-VlMO- 
m-+oo m j^Tj 

Let the weights associated with each link be a(l). 
These weights are a reflection of the relative priority of 
links. We wish to maximize die average throughput (ca- 
pacity) of the network, f(X aV9 )> which is simply the 
weighted sum of long-term average rates: 



m G(T(l),R(l))P(l) 

7W ~ ^ k:k #P(k)G(T(k) y R(l)) +n m 



(1) 



where, nK(i) is the ambient noise power at node R(l). 
Let the raw data rate of link I be X(l). We assume X(l) 
is a linear function of its SIR and the total available band- 
width 1 W. Assuming the Gaussian approximation [6] to 
compute the bit-error-rate (BER), the data rate of a link 
with a tolerable BER of 10~ q and using BPSK modula- 
tion is given by 



X(l) = W'-y(l) where W = 



W 



2glog 2 10 



(2) 



X m (l) = W'U 



G(T(l) 7 R(l))P m {l) 



'£*:**i G{T(k) f R(l))P m (k) + n m 
The long-term average rate of link I is then defined as 

1 m 

x av m) =Bminf — yzx k (i). 



f(X^) = Yla(l)X aV9 {l). 



(3) 



Hie transmitter for each link I is constrained to expend 
an average power of at most P av9 (I) and a peak power of 
atmostP mox (0. 

A Primal Problem 

The objective of this paper is to compute transmission 
schedules for the links in the system so as to maximize the 
system throughput of the wireless network when links are 
subject to peak and average power constraints. The primal 
problem can be succinctly stated as, 



{maxf(X av9 ) subject to P aV9 (I) < P av *(l) 
and 0 < P m (l) < P max } 



(4) 



In fact, we consider a time slotted model, where all slots 
are of equal duration and transmissions begin and end on 
slot boundaries. We let X m (l) and P m (l) be the data rate 
for link I and transmission power for T(l) in slot m. We 
then have 



Note this optimization involves choosing optimal 
power levels for each transmitter in each slot. Further- 
more, the achieved instantaneous throughput in each slot 
is not a convex function of the transmission powers in the 
slot However, using a duality approach below, we reduce 
the problem to a convex optimization problem over a sin- 
gle slot. 

B. Duality Approach 

We exploit a convex duality approach to solve thejari- 
mal problem. Defining a set of dual variables p = 
(0(1), 0(2) . . . P(L a )} for each link, the dual objective 
function is 



g(0) = max{/(X<"*) + f^fi(l)(P avg (l) - P avg (D)} 

subject to 0 <P m <P max (5) 

Note the absence of the average power constraints in 
the definit ion of g(0). For any non-negative vector j8, it is 
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clear tfiafthe maximum value of the objective function in 
the optimization problem (4) is upper bounded by g{0). 
This observation leads to the dual optimization problem, 

mxng{p) subject to {0>O} (6) 

Computation of g{0) involves optimizing over the set 
of powers {0 < Pk < P 1 ™*} in every single slot k. Us- 
ing the linearity of /(.) and the fact that the definition 
of g (0) does not directly incorporate average power con- 
straints, it is clear that g{0) can be evaluated by optimiz- 
ing over a single slot. i.e. g(0) = h(0) where we define 
h(0) below. Defining 



1=1 



(7) 



the dual function over a single slot h(0) is simply 
ft(j3)= max Q{PiJ) 

0<Pi<P max 

Therefore we consider the simplified dual problem 



where a™ = P* v s - P£* b m = Q(P^, 0). We 
have developed an iterative descent algorithm to solve the 
dual problem (8). Exploiting the piecewise linearity of 
h(0), our algorithm converges in a finite number of itera- 
tions. Alternately, by introducing a slack variable, (8) can 
be reduced to a linear program and be efficiently solved 
using interior point methods. Solving the dual problem 
yields 0* and P*«\ such that 

min{h(0) : 0 > 0} = h(0*) - Q{P**J*) for i € [1, JC] 

The vectors P*»* are extreme points of Sp and are 
called the optimal transmission modes. 

IV. Optimal Link Scheduling Policy 
The solution to the dual problem yields us a set of 
optimal transmission modes. The complementary slack- 
ness condition ensures that links with /?*(!) > 0 meet 
their average power constraints with equality. Atleast one 
user will meet his average power constraint with equal- 
ity in the optimal solution. Suppose L eq is defined as 

L eq = 0 : > °}- and let M = \ L .**\* Us- 

ing the principle of strong duality [5], there exists A* = 
(A*(l), A* (2), A*(j£)} such that the equations below are 
satisfied. 



mm 



h(0) subject to {0>O} 



(8) 



III. Solving the Dual Problem 
In order to solve the dual problem, we must first com- 
pute the dual objective function. Consider the polyhedral 
setSp = {Pi : 0 < Pi(0 < P max (lhl = 1,2. ..L G }. 

Let P* 1 * m = 1,2 2 La be the extreme points of Sp., 

Le. each point in Sp can be uniquely represented as a 
convex combination of Pfi** m = 1, 2, . . . 2 L ° . Using the 
feet that Q(Pi , 0) for a fixed 0. is convex in each variable 
Pi (J), the following lemma can be proven 
Lemma 1: 

h{0) = max{Q(P£*,£) : m = 1, 2, . . . 2 L "} (9) 
It is interesting to note that Q{Pu0) is not a convex 
function of Pi, only convex in each coordinate of Pi. 

A Minimizing the Dual Objective function 

The dual objective function g(0) is thepoznftvi.se max- 
imum of 2 La affine functions. Therefore, h(0) is a piece- 
wise linear (convex) continuous function of fi. 



h{0) = max {a m 0 + b m } (10) 



K 

^A^P^'CO = P ot,fl (0 for all links i Gi eg 
k 

J]) A*(i)P**(0 < P aV9 (l) for all links I $ L eq 

i=l 

5^A*(i) = 1 where A*(i)>0 
i=i 

Note that K is not necessarily equal to M + 1, so the 
above system of equations may be under or over deter- 
mined. If K < M + 1, we compute A* by inverting a 
K x K matrix. If K > M + 1, we invert a series of 
Rmax X Rmax matrices, corresponding to different sub- 
sets of optimal transmission modes, where the cardinality 
of the subsets are equal to < M + 1. Note that 

the matrices involved have elements which are either 0 or 
1, which simplifies the computation. We terminate upon 
finding a solution, which is guaranteed by the principle of 

strong duality. 

The optimal policy consists of time sharing between the 
optimal transmission modes, where the fraction of time 
that transmission mode P*' 1 is used is equal to A*(i) for 
each i € [1, jKTJ. In every slot, if a transmitter is included 
in the enabled transmission mode, it transmits at the peak 
power level allowed forjhat transmitter. 
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Simulation results and Discussion 
In this section we compare the performance of our al- 
gorithm with CDMA over the uplink channel for a single 
cell (i.e. single hop) network of 8 nodes. Nodes are scat- 
tered uniformly in a circular region of radius 200 meters. 
The peak power of links are P 1 ^ = 1 Watt Links have 
equal weights i.e. a(l) = |. We compare the throughput 
achieved by the optimal policy against a CDMA policy as 
a function of ambient noise (normalized to the peak re- 
ceived power) and the average to peak power ratio. All 
users in the CDMA policy transmit at their average power 
pavg every s } ot Figure 1 depicts the ratio of the through- 
put of our policy to that of the CDMA policy as a function 
of ambient noise for different values of average to peak 

~- proas pmax , proa* 

power ratio, namely P av * = -7—, —5— and — 5— . 

In low noise regimes, the optimal policy outperforms 
the CDMA policy by at least a factor of 300 for all av- 
erage to peak power ratios considered. The optimal pol- 
icy in low noise regimes schedules good links one at a 
time in TDMA sequence. In this regime, CDMA performs 
poorly due to high mutual interference between transmis- 
sions. As the ambient noise level increases, the optimal 
policy schedules multiple simultaneous transmissions for 
unequal amounts of time. In fact in high noise regimes, 
the optimal policy schedules a large number of simultane- 
ous transmissions; the throughput achieved by the CDMA 
policy therefore approaches that of the optimal policy. 

In low noise regimes, we observe a high improvement 
in throughput with increasing average to peak power ra- 
tios. This is so because the optimal policy schedules good 
links to transmit with increasing frequency (almost equal 
to their peak to average power ratios). At higher ambi- 
ent noise levels, there is relatively little improvement in 
throughput for higher average to peak power ratios. This 
is so because the optimal policy schedules a large number 
of simultaneous transmissions. The optimal policy uses 
lower average power than the specified limit in all noise 
regimes. In fact, in low noise regimes, the average power 
consumed by the optimal policy is significantly less. 
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Fig. 1 . Optimal Policy vs CDMA for low values of ambient noise. 

The set of optimal transmission modes depend on the level 
of ambient noise and average-to-peak power ratios. The 
optimal policy in low noise regimes is to schedule links 
with good channels and low average to peak power ratios, 
in a TDMA sequence. In moderate or high noise regimes, 
the optimal strategy involves scheduling multiple simul- 
taneous transmissions even though they may be in close 
geographic proximity. 
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VI. Conclusion 

We developed an integrated link scheduling and power 
allocation policy for a general multihop network that max- 
imizes network throughput subject to peak and average 
power constraints, assuming constant channels. We be- 
lieve we can extend our algorithms to the case of block 
fading channels. The power control policy is extremal 
and therefore simple. The capacity is achieved by time- 
sharing a small number of transmission modes (< L a J- 1). _ . . 
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Abstract — 

In this paper, we study the problem of joint routing, link 
scheduling and power control to support high data rates 
for broadband wireless multi-hop networks. We first ad- 
dress the problem of finding an optimal link scheduling and 
power control policy that minimizes the total average trans- 
mission power in the wireless multi-hop network, subject to 
given constraints regarding the minimum average data rate 
per link, as well as peak transmission power constraints per 
node. Multi-access signal interference is explicitly modeled. 
We use a duality approach whereby, as a byproduct of find- 
ing the optimal policy, we find the sensitivity of the minimal 
total average power with respect to the average data rate for 
each link. Since the minimal total average power is a con- 
vex function of the required minimum average data rates, 
shortest path algorithms with the link weights set to the link 
sensitivities can be used to guide the search for a globally 
optimum routing. We present a few simple examples that 
show our algorithm can find policies that support data rates 
that are not possible with conventional approaches. More- 
over, we find that optimum allocations do not necessarily 
route traffic over minimum energy paths. 

Index Terms— Wireless Networks, Multi-hop Networks, 
Routing, Scheduling, Power Control. 

I. Introduction and Model 
Broadband wireless networks today are capable of sup- 1 
porting high data rates. Minimizing the total power in 
such systems is of paramount importance not just to in- 
crease its own operational lifetime in the case of battery 
powered devices, but also to coexist symbiotically with 
other systems which share the same frequency spectrum. 
For instance, 802.11 LANs and Bluetooth networks share 
the same unlicensed band and can mutually benefit by lim- 
iting the power of their respective signal transmissions. 

Most currently deployed networks that support high 
data rates use only a "single hop" for wireless commu- 
nication between a mobile user and an access point that is 
connected to the wired infrastructure. However, the use of 
multiple hops to transport data has been shown to enhance 



network capacity [9] [8] [15], and may be necessary due 
to cabling limitations in many environments. 

Consider a network of base-stations interconnected to 
each other through wireless links. Each base-station 
serves as an ingress or egress for the aggregate of traf- 
fic associated with the mobile users in its domain, and 
routes its data through other base stations via multiple 
wireless hops to an access point connected to a wired in- 
frastructure. In this paper, we focus on the transport of 
such aggregated traffic between base stations. We present 
algorithms to compute the optimal routing, link schedul- 
ing and power control policy that consumes minimum to- 
tal average power. While current approaches tend to de- 
couple routing, link scheduling and power control from 
each other, we take an integrated approach. 

Suppose there are N stationary nodes, labeled by the in- 
tegers 1, 2, . . . , N, which correspond to the base stations 
discussed in the previous paragraph. A set £ of Lg = \£ \ 
transmission links, among the possible N{N - 1) links 
between nodes, constitutes a network topology. For now, 
assume that the transmissions on all of these links occur 
concurrently. 

For a given link I = (i, j), the transmitter node s uses a 
signal power P(l). The path "gain" from node i to node 
j is given by G(i 7 j), and models the effects of signal at- 
tenuation due to distance, channel fading and shadowing, 
as well as antenna gain patterns. We assume that the path 
gains G(i 5 j) are constant The transmitting and receiving 
nodes of link I are denoted by T(l) and R(l) respectively. 
The received signal power at node R{1) from the trans- 
mitter T(I) thus given by P(I)G(T(|), *(!))- However, 
signals emanating from other transmitters appear to the 
receiver R(l) as interference, and there is ambient thermal 
noise to contend with as well. The signal to interference 
and noise ratio (SIR ) for link I is defined as 



m G{T{l),R{l))P(l) 

7 ~ IW f WOCT(*).«(0) +*j*i> 



(1) 



where rij is the ambient noise power at node j. 

The capacity of the single wireless channel embodied 
by link I can be modeled as a function of the SIR, 
For instance, the Shannon capacity of a link I over a 
frequency bandwidth W is W\og 2 (l +7(0)» assuming 
Gaussian noise and interference. Throughout this paper, 
we assume that the achieved data rate on link I, X(l), is 
a linear function of the SIR, 7(f). and the total available 
bandwidth W. This is a reasonable approximation to the 
Shannon capacity at small values of SIR, which for exam- 
ple is the case in CDMA networks where all links share 
a common frequency band and have at least a moderate 
processing gain. Our model could easily encompass the 
case where the links may use different frequency bands, 
and in this case the path gains incorporate crosstalk due 
to imperfect bandpass filters. However, for concreteness 
throughout this paper, we shall assume all links share the 
same frequency band of width W. Assuming the Gaussian 
approximation [12] to compute the bit-error-rate (BER), 
the data rate of a link I with a tolerable BER of 10~* and 
using BPSK modulation is given by 



X{l) = W'~y{l) where W = 



W 



2$log 2 10* 



(2) 



. Suppose that the desired data rate on each link I is C [l) 
bits per second, i.e* 

X{l)>C(l).lGL E . (3) 
Using (1) and (2), we can rewrite (3) in matrix form as 

P>PP+b, (4) 
where P = (P(l), .... P{L £ )) T , 

K ' G(R(l),T{l))W" 

b = (6(l),...,6(X,f)) T , and P is the L £ x L E matrix 
with entries 



f ^ 

*■</,*) = { 0 < 



GWQ.Tffl)C(l) 



,iffe^/ 
,iffe = J. 



If F is a stable matrix, Le. all eigenvalues of F are 
strictly inside the unit circle, then (I — F)"" 1 exists and 
has non-negative elements. Thus, in this case, we can re- 
arrange (4) and premultiply the resulting matrix inequality 
by (I- F)"* 1 to obtain . 



<£P --TP »U» <U& t "IS » W w > JL? £jf . 



P > (I-F^b 



(5) 



Hence, we can view (I — F)~ x b as the minimal power 
vector that supports the network topology defined by links 
£ and their associated required data rates. Foschini and 
Miljanic [7] proposed a simple distributed algorithm to 
compute the minimum power vector in a cellular network 
context Subsequently, Bambos, Chen, and Pottie [1] con- 
sidered a distributed link admission control algorithm for 
the case of multi-hop networks, in the framework of [7], 

If the eigenvalues of F are on or outside the unit circle, 
there does not exist any power vector P which supports 
the required data rates. If the eigenvalues of F are in- 
side the unit circle but close to the boundary, the minimal 
power vector (I — F)~ x b will be very large. Even if the 
eigenvalues of F are well inside the unit circle, using the 
minimal power vector indicated in (5) may be inefficient 
This is illustrated in the following example. 

Example: Consider a square network of 4 nodes and 
2 links {(1,2), (3,4)}. Suppose the relevant path gain 
values ares G(l,2) = C?(3,4) = G 0 and £(1,4) = 
C?(3,2) = Go/2, corresponding to omnidirectional an- 
tennas and inverse square law propagation loss. Let the 
required data rates be X({1, 2}) = X({3, 4}) = and 
let the ambient noise no be the same at both the receivers, 
specifically no = Gq. In this case,' it is readily verified 
that (I - F)~ x .b = (0.67, 0.67) T , i.e. each transmitter is 
required to continuously transmit at a power of 0.67 Watt. 
On the other hand, if instead the transmitters take turns 
transmitting at a power of 1 Watt, such that only one trans- 
mitter is active at any time, the peak data rate that can be 
achieved is twice that of the data rate achieved previously, 
due to a lack of interference. Thus, if the transmitters are 
active for 50% of the lime, the same average data rate of 
^ is achieved, but this consumes an average power of 
only 0.5 Watt, which is 33% less than before. 

This example suggests that it can be inefficient to op- 
erate with all links in a wireless network active concur- 
rently. Indeed, in most wireless networks, links operate 
in a half-duplex mode to limit interference. We can es- 
sentially eliminate interference by scheduling only links 
which are very far apart, as is proposed for example in 
[61[10]. However, in order to achieve high data rates, this 
may require unacceptably high peak transmission power. 
Constraints on peak transmission power may necessitate 
the use of concurrent transmissions in this case. 

Recently, Elbatt and Ephremides [5] proposed a 
scheduling and power control strategy that selects candi- 
date subsets of concurrently active links, and applies the 
_ power control algorithm in [7] in order to find the mini- 
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mal^r vector (I - P) -1 b for each candidate subset 
The objective of their optimization problem is to maxi- 
mize network throughput by supporting as many links as 
possible with a given fixed SIR, while remaining energy 
efficient. They propose a heuristic to identify large subsets 
of concurrently active links that can be efficiently sup- 
ported, and evaluate the throughput of the heuristic via 
simulation for a uniform traffic demand on links. 

In this paper, we propose a systematic optimization for 
a problem formulation that is motivated by similar ob- 
jectives. In particular, we seek to find subsets of simul- 
taneously active links, herein called transmission modes, 
as well as the associated transmission powers, in order to 
minimize the total average transmission power expended 
across the network. The constraints of this optimization 
problem are that each link in the network must have an 
average data rate no less than a given prescribed value, 
and a peak transmission power for each node. This prob- 
lem statement is formalized in Section n, and a duality 
approach for finding the optimal scheduling and power 
control policy is presented. 

The required minimum average data rates on each link 
are ultimately determined by the average rates at which 
traffic is generated by users and the routing algorithm that 
is used. In wireless networks it is possible to reconfigure 
the data rates of the links on a fast time scale in response 
to changing traffic and channel conditions. In this paper, 
we shall also consider the problem of routing, and hence 
determination of the required data rates on each link, for a 
given traffic demand rate matrix. We shall address this in 
Section m, exploiting the duality approach explained in 
Section II. 

After presenting the details of our approach, we discuss 
several illustrative examples in Section IV. Qualitatively, 
we find that in high ambient noise regimes, a large number 
of links in close geographic proximity are active simulta- 
neously in optimal policies. In addition, we find that route 
diversity allows nodes to split traffic over multiple paths 
and support higher throughputs than could otherwise be 
achieved. This is somewhat surprising since all links share 
a common bandwidth. However, pushing the data rates to 
the maximum possible limit comes at price: As the data 
rates increase, the network tends to become less energy 
efficient 

We conclude the paper with some closing remarks in 
Section V. 

n. Optimal Scheduling and Power Control 

In this section, we formally specify the optimization 
problem discussed in the previous, section, and present a 



duality approach for solving it We then discuss an ap- 
proach to limit computational complexity. 

For simplicity of exposition, we divide time into slots, 
each of equal duration and indexed by the positive inte- 
gers. Transmissions begin and end on slot boundaries. 
Generalizing the notation introduced earlier, let -X* m (i) 
and P m (J) be the data rate for link J in slot m, and trans- 
mission power for the transmitter T(i) for link I in slot ro, 
respectively. Let P m - (P m (l),Pm(2), - . - , Pm{Ls)) be 
the network power vector for slot m. Let P max (i) be the 
maximum transmission power for node i. Let S (i) be the 
links in £ that originate at node «. 

Each node must conform to the peak transmission 
power constraint in every slot: 

0 < PmV) < P max {i) and (6) 

0 < P m (0 ,forallTO>l. 
Using (2) and (1), the achieved data rate for link I in 
slot m is 

v m n-r g(r(0,fl(0)Pm(0 , 

XmW - W ^ kzk #G{T{k)Ml))Pmik) + 
The long-term average rate of link I is then defined as 

1 m 

X av9 {l) = liminf — Y]X k (l). 

k=l 

For each link /, let C(l) be a given minimum acceptable 
average data rate, i.e. we must have 

X avg (l) > C(l) , for all I € 5. (7) 

Define the required minimum average rate vector as C = 
(C(1),C(2),...,C^)). 
The average power consumed by the transmitter tor link 

l,P aV9 (!),is 

1 m 

P™ 9 (1) = limsup — rPfc(l)- 
m->oo rn f— ' 



Define the average network power vector as P aV9 = 
(P<"tf (1),P<"*(2), _ , P*»9(L S ))- 

There may or may not exist a sequence of network 
power vectors P u ft, - - - that satisfy (6) and (7). If there 
does exist a sequence of such network power vectors, our 
aim is to minimize a linear function of P aV9 - An example 
of such a linear function is simply the total average power. 
We will consider the slightly more general cost function 



M^) = 2>(0P°*'(0, 



(8) 



whereof) is a positive weight For example, hiP™ 9 ) 
might equal the average power incident at a given point in 
space, or be equal to the total average power emanating 
from a closed boundary surrounding the network. Alter- 
natively, a(l) might model the cost of a battery that drives 
transmitter T(J). 

A. Primal Problem 

One of the main aims of this paper is to present an 
approach Jo compute a schedule of network power vec- 
tors P\,P2, • - . • satisfing the peak power constraints (6) 
and the minimum average rate constraints (7), such that 
h(P aV9 ) is minimized. This primal problem is succinctly 
stated as 



mm{fc(P ot *)} subject to (6) and (7). 



(9) 



Note this optimization involves choosing optimal 
power levels in each slot for each transmitter. Using a du- 
ality approach below, we reduce the problem to a convex 
optimization problem over a single slot 

£. Duality Appwacli 

Let the value of the optimal cost in the problem (9) as a 
function of Q be denoted by H(C). If for a given value of 
C, no schedule of network power vectors exists satisfying 
(6) and (7), define H(C) = +00. 

Note that H{C) is a convex function of the vector &. 
To see this, suppose that <? = AC?i 4- (1 - A)C7 2 , where 
0 < A < 1. Let 7Ci be an optimal policy that supports 
the average date rates specified by i.e it achieves mini- 
mum cost H(Ci). Roughly speaking, we can create a new 
policy 7r by 'time sharing" between the policies tti and 
7T2, such that policy 7ri is used a fraction A of the time and 
policy 7r 2 is used a fraction 1 — A of the time. The pol- 
icy TV supports average data rate rates that are at least as 
large as C. Furthermore, by the linearity of fr(-), the cost 
of policy 7t is \H(Ci) + (1 - X)H{C 2 Y By definition 
of H{C), the cost of policy 7r is at least H(C). Thus we 
have H(C) < XH(Ci) + (1 - X)H(C 2 ). 

Define a set of dual variables 0 = 
{/3(1),)9(2),... ,0(L e )} for each link. Define a 
•^potential" function V according to 

v(P,d y x) = h(P) + ][>(o [c(i) - x(i)} 

lee 

The dual objective function is defined as 



g0) =mln{V(P av9 7 C } X av9 )} subject to _(6) 



*L2l «US» "IT -3 jt -y» W H«i» ^ ^3 



Note the absence of^the minimum rate constraints (7) 
in the definition of g{0). For any non-negative vector 0 9 
using (7), note that 

H(C) > Inm{^(P* v ^C,J? a ^)},subjectto(6)anda) 
> mhx{V(P av9 , C,£™ 9 )} , subject to (6) 

= 90) * (10) 

Thus, for any non-negative vector 0 y the minimum 
value of the objective function in the optimization prob- 
lem (9) is lower bounded by 90). This observation leads 
to the dual optimization problem, 



(ID 



A geometric interpretation of the inequality (10) can be 
obtained by considering the surface of H (•) in (Ls + 1)- 
dimensional Euclidean space. In particular; H(S) is 
lower bounded by considering a supporting hyperplane 
that is normal to the vector (j8, 1). Since H(*) is convex, 
there exists a supporting hyperplane which passes through 
the point (C,H(C)). This is simply the classical argu- 
ment used to prove that there is no "duality gap," [2] i.e. 



H(C) = maxteO?) :0>8] 
= 9(0*)* 



(12) 



Here, the supporting hyperplane passing through the point 
((?, H [6)) is normal to the vector (jS*, 1). From this ge- 
ometric interpretation, it can be seen that the components 
of the optimal dual variable vector 0* represent a sensitiv- 
ity of the optimal cost with respect to a perturbation in the 
minimum average data rate for a link, i.e. 



r (0 = 



3H{C) 
dC(l) 



(13) 



Computation of g[fi) involves optimizing over all 
schedules of network power vectors satisfying the peak 
power constraint in every slot However, since the poten- 
tial function V(P,d,X) is linear in P,C, and X, it fol- 
lows that g 0) can be computed by an optimization over 
a single slot, i.e. 

90) = min{K(P, C,X):Pe S P } , 

where X = X(P) is defined as in (2)juid (1), and Sp is 
defined as the polyhedral set Sp = {P : Y^ieeu) F ( l ) < 
P^WforalltandO < P(l) foralU € £ }. Let M 
be the number of extreme points of Sp, and let P^f*, 
m = 1,2,... Af denote the extreme points of Sp. By def- 
^ch^oint in Sp can be represented as a convex 
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combination of P^ 1 , m = 1, 2, . . . M. An upper bound 
on M is 2 Le . If Mi is the number of links emanating from 
node i> then M = nj^l + Mi). 

Even though it can be verified thatV^P, C, X (P)) is 
not a concave function of the vector P, it is easy to show 
by taking derivatives that V(P 1 C 7 X(P)) is a concave 
function of P(i) for each I € 5. A standard convexity 
argument then shows that 

9(0) - minO^S <?, j£(P£*)) : 1 < m < M} . 

Thus, the complexity of computing g(0) is 0(M). In 
fact^note that V(P^ t y C y X(P^ t )) is an affine function 
of 0, which implies that the dual objective function g(0) 
is the pointwise minimum of M affine functions, Le. 

g(0) = min {a' m /9 + 6 m > . (14) 

We have developed an iterative ascent algorithm to 
solve the dual problem (1 1). Exploiting the piecewise lin- 
earity of g(0), our algorithm converges in a finite number 
of iterations. Alternately, by introducing a slack variable, 
the dual problem (1 1) can be reduced to a linear program 
and be efficiently solved using interior point methods. 

If the required minimum rate vector C is infeasible, i.e. 
H(C) = -t-oo, it can be seen that g(0) is unbounded from 
above over the set of all (3 > 0. If C is feasible, it is easy 
to see that 

N 

H{6) <^P ma *(i). 

i=l 

Thus we can safely declare C infeasible if we find a value 
of 0 such that g{$) > J^L X P™*®. 

G Computing the Optimal Policy 

Solving the dual problem (11) yields the optimal dual 
variable vector /S* and extremal network power vectors 
P*>\ such that 

g{0*) = mi*{g(0):0>O} 

= V(P^\C,X(P^)) . (15) 

The vectors P* ji are extreme points of Sp that repre- 
sent optimal transmission modes and the power levels of 
the transmitters in those transmission modes. In other 
words, assuming H(C) is finite, an optimal schedule of 
network power vectors exists that consists solely of these 
extremal network power vectors P**\ This means that 
in the optimal policies we find, in each slot, every node 
is either not transmitting at all, or is transmitting at the 
maximum possible peak power to exactly one receiver. 

i 
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Moreover, since a hyperplane in (Le + 1) dimensional 
Euclidean space is determined by Le + 1 linearly inde- 
pendent points contained within it, an optimal policy can 
be constructed that consists of at most Le + 1 extremal 
network power vectors P* fl . Let K be the number of ex- 
tremal network power vectors P*»* such that (15) holds. 
Define K = mii^lC, Le + 1}- We thus have 

H(C) = g(0') = V(P^ 7 C y X(P^)) , 

for all s such that 1 < i < K. 

It is easy to see that in any optimal policy (assuming 
H(G) is finite), the average rate on each link f, X av ^(l) y 
is exactly equal to C (I) . Otherwise, we could decrease the 
total average power further without violating any rate con- 
straints thus contradicting the optimality of the solution. 
Let X*>*(1) denote the rate of link I corresponding the 
optimal extremal network power vector P*»\ Assuming 
H(C) is finite, it can be shown that (12) guarantees the ex- 
istence of a "weight vector" A* = {A*(l), A* (2), A*(£)} 
such that the equations (16) and (17) below are satisfied: 

k 

]T A* (i)X** (I) = C(l) for all links I € 6 (16) 

k 

53 A* (i) = 1 where A*(t) > 0 (17) 
»=i 

Note that K < Le + 1, so the above system of equa- 
tions is always overdetermined. There must be K + 1 of 
these equations which are linearly independent. Once we 
have identified such a set of equations, we solve them by 
inverting a (K 4- 1) x (K + 1) matrix. 

The value of X*(i) indicates the relative frequency at 
which the extremal network power vector P*»* is utilized 
in an optimal policy. This then determines a family of op- 
timal policies, which can be constructed by alternating be- 
tween the extremal network power vectors in accordance 
with the weight vector A*. A specific optimal policy, for 
example, can be constructed using a psuedo-random num- 
ber generator to select which extremal network power vec- 
tor is used in each slot, such that the network power vector 
P*» i is used with probability A*(i). Alternatively, a spe- 
cific optimal schedule can be constructed in such a way as 
to provide "smooth service" to each link, though this topic 
is beyond the scope of the present paper. 

Z>. Reducing Complexity 

The algorithm developed to solve the dual problem con- 
siders all M transmission- modes. However, in the opti- 
mal schedule we have no more than Le + 1 transmission 



modes^uice M can be as large as 2 Le , the complexity 
of minimizing the dual objective function can be expo- 
nential. Thus, reducing the set of possible transmission 
modes considered can greatly reduce complexity. A natu- 
ral consequence of our algorithm for a network of nodes 
with omni-directional antennas is that, no node transmits 
or receives data at the same time in the optimal sched- 
ule due to high self interference. This observation can 
reduce the number of candidate transmission modes sig- 
nificantly. In general, in many cases it is possible using 
heuristics to eliminate many "obviously inefficient" trans- 
mission nodes from consideration, since they are unlikely 
to be used in an optimal policy. This idea is similar in 
spirit to the approach used in [5] to construct a set of 
**valid" links. The resulting policies we find can be shown 
to be optimal among all policies which are constrained to 
use only the set of chosen candidate transmission modes. 

III. ROUTING AND CAPACITY ALLOCATION 

In the previous section, we presented an algorithm that 
finds an optimal link scheduling and power control policy 
to support a given traffic rate on each link in the network. 
In this section, we build on this to construct a policy which 
allocates traffic rates on each link in accordance with a 
routing algorithm, given a traffic matrix which describes 
the rate at which traffic is to be carried between specific 
source-destination pairs. The routing algorithm is guided 
by the computation of the minimum cost associated with 
optimal scheduling and power control for supporting traf- 
fic rates on each link that have already been allocated, as 
described in the previous section. 

For example, consider a series of iterations, where at 
a given iteration we wish to support a small amount of 
additional traffic from one of the given source-destination 
pairs. Initially, no traffic is supported, i.e. the required 
minimum data rate on each link is zero. We compute the 
optimal values for the dual variables, 0*. Using (13), we 
can estimate the cost of supporting additional traffic on 
each link using a first order Taylor approximation. For 
example, suppose we wish evaluate a given route r. Let 
v be the L& dimensional vector, such that the component 
r (I) = 1 if link I is included in the route f, and r(2) = 0 
otherwise. If we route an additional e units of traffic along 
route r, the required data rate vector will increase from C 
to C+er. Using (13), the cost of supporting the additional 
e units of traffic along route r can be determined from the 
approximation 

H{C + er) = H{C) + e £ 0*<|)r (£) . 

i 
I 
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Therefore, in order to determine the best possible path 
along which an additional increment of traffic should be 
routed, we can use a shortest path algorithm (such as Di- 
jkstra's algorithm or the Bellman-Ford algorithm), where 
we set the weight for link I to the sensitivity value 0*(O- 
After we have identified the shortest path, we attempt to 
support the additional traffic along this path. If the new 
value of C is infeasible, we either reduce the step-size e 
until the additional traffic can be supported, try support- 
ing additional traffic for another source-destination pair, 
or otherwise reach a prespecified termination condition. 

In the examples, we will present below, it is apparent 
that jHT(-) is piecewise linear, which results in "sharp cor- 
ners". When we use the above algorithm, it is possible 
that the link sensitivities can change abruptly. This can be 
partially addressed by choosing the step size e sufficiently 
small. We can also '"backtrack" by re-adjusting the rate 
allocations after each iteration, if necessary, to shift traffic 
towards shorter paths. 

Indeed, since H (•) is convex, we can use the optimal 
routing framework developed by Bertsekas et al. (see [3]). 
In this class of algorithms, any feasible routing is first cho- 
sen. By incrementally adjusting the routing in an iterative 
manner, flows are shifted towards shortest paths at each it- 
eration, and eventually a globally optimum routing is ob- 
tained. There are many possible variations, for example 
regarding whether flow is shifted simultaneously for ev- 
ery source destination pair, or whether flow is shifted one 
at a time for each source destination pair. In the present 
context, we are faced with a more difficult problem, since 
apparently J?Q is not differentiate everywhere. In addi- • 
tion, since evaluating the function £f(-) is not trivial, it is 
more difficult to ascertain whether a given routing is fea- 
sible, i.e. the capacity of a given link depends on the rate 
carried by other links. Nevertheless, our initial numerical 
experiments that are based on (13) and framework devel- 
oped by Bertsekas et. al. suggest that such an approach 
may be quite useful in practice. 

As expected, since H (») is convex, we find that fi*(l) 
increases after we allocate additional traffic on link L In- 
terestingly, we find that sometimes will decrease 
when we allocate additional traffic on a link V ^ /. This 
corresponds to negative cross partial derivatives of H (•), 
and does not contradict the fact that Jf(-) is convex and 
non-decreasing. As a result of this, routes which initially 
appear unattractive sometimes become more attractive af- 
ter traffic is added on other links. As we will see, this 
phenomenon can lead to traffic from the same source- 
destination pair being routed over multiple paths, which 
is somewhat counterintuitive since all links share a com- 
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\?l^^nark that the approach of using optimal values 
of the dual variables as link weights in shortest path rout- 
ing has also been used before in the context of revenue 
maximization in wired networks, by Mitra et. al- [11]. 
This work focuses on the asymptotic behavior of the link 
sensitivities with respect to uniform increases in demand 
rates. The authors of this previous work also observed the 
optimality of route diversification for price elastic traffic 
associated with maximizing network revenue. 

The potential energy efficiency of a given link can be 
evaluated, assuming that all other links are idle. For our 
model, the minimum transmission energy required per bit 
on a link (t, j) is inversely proportional to G{i 7 i). A min- 
imum energy path is defined as a path which has the small- 
est possible energy required to transport a fixed amount of 
data, irregardless of the delay associated with transporting 
that data. Previous work on routing in wireless networks 
(for example see [14][13]) is focused on energy efficiency 
without regard to delay, i.e. such work proposes using 
minimum energy paths. 

In our framework, we find that the paths over which 
traffic is allocated in some optimum routings does not al- 
ways correspond to using minimum energy paths. How- 
ever, when required data rates are sufficiently low, our 
routing framework will apparently always use minimum 
energy paths. This is illustrated by example in the next 
section. 

IV. Resource Allocation Examples 

In this section, we illustrate the gains of using our wire- 
less resource allocation framework through a series of net- 
work topologies. We focus on simple examples in an at- 
tempt to gain insight from optimal allocations. 

A. Example 1: String Topology 

A string topology consists of a row of nodes connected 
by means of directed links as shown in Figure 1. This 
topology represents a "transport" network carrying data 
over long distances. Node 1 is the source of data and node 
5 is the sink. Each node in this example has an omnidirec- 
tional antenna, and all nodes have identical peak transmis- 
sion power constraints of 1 Watt The path-loss between 
nodes was modeled as an inverse square of the distance 
for all the examples in this paper, i.e. G(i, j) = g{Ap - 
The ambient noise power at all the nodes is assumed to 
be constant We compare the maximum throughput of the 
single session between node 1 to node 5 achieved by our 
link scheduling policy to a TDMA scheduling policy. The 
TDMA policy schedules exactly one link at a time, ac- 
tivating each link at maximum power for ^fourth of the 



time. We plot the ratio of throughputs achieved by the 
our link scheduling algorithm to TDMA with increasing 
values of ambient noise in Figure 2. The ambient noise 
is normalized with respect to the peak received power 
(from a single transmitter) in the network and is scaled 
by log 2 (.) (for all graphs in this section). The maximum 
throughput achieved by our policy is higher than TDMA 
for all values of ambient noise. In low noise regimes, 
TDMA is a near-optimal policy. However, as the level of 
ambient noise at nodes increases, the gains in throughput 
of our policy increase significantly. The optimal policy 
schedules multiple transmissions at the same time, even 
in fairly low ambient noise regimes, thereby exploiting 
spatial reuse and thus outperforming TDMA. This effect 
is more pronounced in moderate and high ambient noise 
regimes where the optimal policy schedules concurrent 
transmissions even though they are in close geographic 
proximity. 

To get a sense of how energy efficient the optimal pol- 
icy is as it supports increasing traffic demands, we plot 
the total average power with increasing traffic rates for 
the session between node 1 and node 5 in Figure 3. For 
this experiment, we fix the level of the logarithm of the 
normalized value of ambient noise to -0.67. This value 
corresponds to the 3 rd data point on the x-axis in Figure 
2. WealsosetW = 10 7 . 

Figure 3 indicates that the total minimum average 
power is a linear function of the traffic demand for av- 
erage rates below 4 Mbits/sec. Below 4 Mbits/sec, our 
scheduling policy reduces to TDMA, with only one link 
transmitter active at any time. As the traffic demand of 
the session increases beyond 4 Mbits/sec, the total average 
power increases at a faster rate. For traffic load between 
4 Mbits/sec and 4.98 Mbits/sec, the transmission modes 
in the optimal policy are given by [{(1,2)}, {(2,3)}, 
{(3, 4)} and {(1, 2), (4, 5)}]. For traffic loads above 4.98 
Mbits/sec to the maximum throughput possible, the opti- 
mal transmission modes are [{(1,2)}, {(2,3)}, {(3,4)}, 
{(1,2), (3,4)} and {(1,2), (4,5)}]. We see that in order 
to support high traffic loads, the optimal link scheduling 
policy activates a large number of links simultaneously. 
Note that the transmission modes in the optimal policy 
are half-duplex. Therefore, transmissions in the optimal 
schedule cannot consume a total average power of more 
than 2000 mW for this topology. 



B. Example 2: Diamond Topology 



In this experiment, we- consider an asymmetric dia- 
mond topology as_ shown in Figure 4 to illustrate the sub- 
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Fig. 2. String Topology: RaUo of maximum throughput for optimum 
scheduling VsTDMA. 




~~ 1 2 3 * s 

Increasing Traffic toad (Mbits/sec) 




Rg. 3. String Topology: Total Minimum Average Power in milliwatts 
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Fig. 5. Diamond Topology: Ratio of throughputs achieved by using 
multiple paths to using a single path 1 2 -> 4. 

stantial increase in throughput by splitting traffic over 
multiple routes including paths that are energy inefficient 
Nodes in this topology are equipped with omnidirectional 
antennas, and the peak transmission power each node is 
fixed at 1 Watt. Node 1 is the only source of data and 
node 4 is the sink. We assume G(l, 4) = the path 

loss between all other nodes is given by the inverse square 
law of distance. We consider our integrated routing and 
scheduling algorithm over this diamond topology, where 
all possible routes from the source to the destination are 
allowed except for the single hop route (1,4). For com- 
parison purposes, we also considered our scheduling algo- 
rithm where only the links (1, 2), and (2, 4) are allowed to 
be used, which corresponds to routing over a single mini- 
mum energy path 1 2 4. 

We compare the ratio of throughputs for these polices 
in Figure 5. Clearly, using multiple paths yields higher 
throughputs for all values of ambient noise. The increase 
in throughput by splitting traffic over multiple paths is sig- 
nificant even for moderate levels of ambient noise. This 
result is somewhat surprising, since all links are shar- 
ing a common bandwidth and we have a per-node peak 
power constraint. An intuitive explanation of this is that 
by splitting the traffic over both paths, the transmission 
modes {(1,2),(3,4)} and {(2,4), (1,3)}, which do not 
have much interference, can be alternated in time. This 
results in substantial delivery of data to the destination in 
every slot If we only allow the path 1 -» 2 4, in or- 
der to avoid self interference, the links (1,2) and (2,4) 
must be active one at a time, so that node 4 will only be 
delivered data in every other slot 



To illustrate energy efficiency issues, we plot the total 
minimum average power with increasing traffic demand at 
node 1 in Fieure 6. The logarithm of the normalized value 
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Rg. 6. Diamond Tbpology: Total Minimum Average Power in milli- 
watts Vs Increasing traffic demands in Mbits/sec 

of ambient noise for this point is 0.45 and corresponds 
to the 3 rd data point on the x-axis in Figure 5. We set 
W' = 10 7 . The maximum data rate achieved by TDMA 
using the minimum energy path 1 -> 2 -> 4 is 3-65 
Mbit/sec. In comparison, our policy is capable of support- 
ing a data rate of up to 5.05 Mbits/sec by using bom the 
minimum energy path as well as the non-minimum energy 
path. For traffic loads below 3.65 Mbits/sec, our policy is 
identical to TDMA over the path 1 -+ 2 -» 4. As the traf- 
fic load increases beyond 3.65 Mbits/sec, the optimal pol- 
icy uses starts using both paths by scheduling transmission 
modes {(1,2), (3,4)} and {(1, 3), (2,4)} in addition to 
TDMA transmission modes. The average power required 
to achieve higher traffic demands beyond 3.65 Mbits/sec 
increases with a greater slope as can be seen in Figure 6. 
The maximum throughput of 5.05 Mbits/sec is achieved 
by splitting traffic by using both paths and scheduling 
transmission modes {(1, 2)}, {(2, 4)}, {(1, 2), (3, 4)} and 
{(1,3), (2, 4)}. 



C Example 3: Hierarchical Topology 

A hierarchical topology as shown in Figure 7 represents 
a complex topology such as an access network, where the 
number of nodes in each tier decrease as we get closer 
to the access point node 6. This experiment highlights 
the significant gains in throughput by jointly routing and 
scheduling to an alternative policy that optimally sched- 
ules links over pre-determined minimum energy routes. 
Nodes 1, 2 and 3 are sources of data and node 6 is^the 
sink. We assume G(l, 6) = ap^r and G(2, 6) = a^jr; 
the path loss between all other nodes is given by the in- 
verse square law of distance. The minimum energy paths 
for the source nodes are: {1 -» 3 -¥ 6}, {2 -> 3 -*• 6} 
and {3 -)• 6} respectively. Our integrated routing and 
scheduling algorithm is allowed to utilize all possible 
paths to route their data to node 6. We consjrainffie .rates 



Rg. 7. Hierarchical Topology: Access Network 

of each source to be identical. We plot the ratio of the 
total maximum throughput of our routing policy to that 
achieved by using only the minimum energy paths, as a 
function of ambient noise in Figure 8. 

We find that at low noise levels, our policy supports a 
throughput at least 30% more than the data rate achieved 
using minimum energy paths alone. As the level of am- 
bient noise increases, our policy supports a significantly 
higher throughput. The base policy schedules links in a 
TDMA sequence for all values of ambient noise. The 
optimal policy schedules concurrent links even for low 
levels of ambient noise. In the optimal policy, nodes 1 
and 2 route their traffic over non-minimum energy paths 
1 4 6 and 2 -» 5 -> 6 respectively, apart from the 
minimum energy paths, for all values of ambient noise. 
For low ambient noise, the optimal policy schedules the 
transmission modes {(1,3)}, {(2,3)} and {(3,6)} ac- 
cording to TDMA but also allows modes {(1,4), (5,6)} 
and {(2, 5), (4, 6)} to be active for reasonable fractions of 
time. As the level of ambient noise increases, the opti- 
mal policy splits node 3's traffic over paths 3 -» 4 ->• 6 
and 3 5 -i- 6. In fact, for sufficiently high ambi- 
ent noise, the optimal policy schedules a large number 
of concurrent links, the dominant transmission modes be- 
ing {(1,4), (5, 6), (2, 3)} and {(2,5), (4,6), (1,3)} and 
{(1,4), (2, 5), (3, 6)}. 



V. Conclusion 
We have developed an integrated routing, link schedul- 
ing and power allocation policy for a general multihop 
network that minimizes the total average power consump- 
tion to support minimum average rate requirements per 
link. Our policy can support higher throughputs than with 
conventional approaches to radio resource allocation, at 
the expense of decreased energy efficiency. Our policy re- 
quires time synchronization between transmitters, and re- 
quires that channel conditions remain constant over sev- 
eral time.slots. It appears to be most appropriate in the 
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Fig. 8. Hierarchical Topology: Ratio of maximum throughputs 
achieved by using multiple paths to that using minimum energy paths 
only 



context of providing a means for wireless interconnection 
of fixed stationary nodes, such as wireless access nodes. 

We find that our optimal power control policies are such 
that each node is either transmitting at peak power to a sin- 
gle receiver or not transmitting at all, which makes it sim- 
pler to implement This also means that nodes can go to 
"sleep" much of the time, when they are not transmitting 
or receiving, which potentially conserves a great deal of 
energy [4]. The optimal link schedule time-shares a small 
number of optimal subsets of links (< Ls + 1) in order 
to achieve the required data rates. The collection of op- 
timal transmission modes depend on the level of ambient 
noise and the capacity allocations on links. The optimal 
policy for low required data rates or in low ambient noise 
regimes schedules links in a TDMA sequence. In moder- 
ate or high noise regimes, or to achieve higher data rates, 
the optimal strategy involves scheduling multiple simul- 
taneous transmissions even though they may be in close 
geographic proximity. 

A byproduct of solving the scheduling and power con- 
trol problem yields link sensitivities that tell us the change 
in the minimal total average power with respect to a per- 
turbation of the required data rate on a link. The mini- 
mal required average power is a convex function of the 
required data rates on each link. This leads us to a rout- 
ing algorithm that is based on previously developed non- 
linear programming techniques. 

We find that using minimum-energy paths to route traf- 
fic and TDMA to schedule links is near-optimal in low 
noise regimes and in the case where the required aver- 
age data rates are small. As the level of ambient noise 
increases, or the required average data rates increase, our 
optimal routing policy exploits the use of multiple paths to 
support higher data rates, and the optimal scheduling pol- 
icy makes use of concurrent transmissions on links. We 
find that non-minimum energy paths can be exploited to 



increase throughput, despite the fact that all links share a 
common bandiwdth. 

We remark that it is easy to generalize our problem for- 
mulation in a number of directions. For example, we can 
specify a maximum peak power constraint for an arbitrary 
collection of nodes, rather than per each node. This may 
be useful to limit signal emissions from an entire network 
or subnetwork. 

One topic we are currently investigating is the use of 
hierarchical scheduling, whereby links are partitioned into 
clusters that are geographically remote. The advantage of 
this is that the clusters can be scheduled independently, 
with only a loose coupling between different clusters. The 
signal emissions from a given cluster can be modeled as 
"ambient noise" for purposes of scheduling transmissions 
in another remote cluster. 
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